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One of the major frontiers of lattice field theory is the inclusion of light fermions in simulations, particularly 
in pursuit of accurate, first principles predictions from lattice QCD. With dedicated Teraflops-scale computers 
currently simulating QCD, another step towards precision full QCD simulations is underway. In addition to 
ongoing staggered and Wilson fermion simulations, first results from full QCD with domain wall fermions are 
available. After some discussion of work toward better algorithms, simulations completed to date will be discussed. 



1. INTRODUCTION 

A major challenge in realizing the full potential 
of the non-perturbative regularization provided 
by lattice field theory is the presence of fermion 
fields in simulations. In many cases, such as QCD 
at finite density and the Hubbard model, the 
fermions require an improvement in algorithms 
before large systems can be simulated. For QCD 
at zero chemical potential, simulations with dy- 
namical fermions will continue to become bet- 
ter controlled by faster computers, even without 
(hoped for) theoretical improvements. 

There is over a decade of experience with full 
QCD simulations, although only recently has 
enough data become available for some extrap- 
olations to the continuum As with quenched 
simulations questions of what are satisfactory vol- 
umes and lattice spacings, what are reliable sim- 
ulation lengths and what are the best extrapola- 
tions need to be answered. Removing the uncon- 
trolled truncation of the quenched approximation 
is vital to achieving precision results, but answers 
to the above are very important. 

After a brief overview of the work on algorithms 
presented at this conference, I will focus on re- 
sults for QCD, particularly the low lying hadrons. 
There is evidence the quenched hadron spectrum, 
in the continuum and chiral limits, differs from 
nature 11 . For fuU QCD, the light hadron spec- 
trum is interesting in its own right, provides the 
basis for a determination of light quark masses 
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and is a useful testing ground for exploring fi- 
nite volume effects and simulation lengths. Also, 
hadronic states differing only by their parity are 
degenerate unless chiral symmetry is broken, so 
studying these other parity states in the spectrum 
may show sensitivity to the inclusion of the quark 
determinant. 

Most of the techniques developed in the 
quenched approximation for measuring other ob- 
servables can be easily adapted to full QCD sim- 
ulations. Many groups are reporting on effects 
of dynamical fermions outside of the light hadron 
spectrum {/b, etc.). Please see the reviews of 
these areas for information. 

2. ALGORITHMS 

Three groups reported work on the algorithms 
used to produce the Markov chain in numerical 
simulations. 

2.1. Meron Cluster Algorithm 

Wiese and collaborators (MIT) have tack- 
led the long standing problem of algorithms for 
fermion systems where the weights in the path 
integral are complex. Computational times expo- 
nential in the volume are needed if the weight is 
made real by taking its absolute value and mov- 
ing the phase to an observable. By working with 
a continuum Euclidean time formulation, where 
fermions are represented by their world lines, 
their algorithms form clusters whose flip gives a 
definite sign. A second improvement insures that 
clusters which cancel through a sign flip and those 
that don't are generated with similar probabili- 
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ties. This algorithm, with its attendant improved 
estimator, has allowed them to simulate various 
fermionic systems that are intractable with con- 
ventional approaches. Application of these ideas 
to coupled gauge-fermion systems is being inves- 
tigated. 

2.2. Truncated Determinant Approach 

Duncan, et. al. Q| have been studying an al- 
gorithm where the fermion determinant is split, 
in a gauge-invariant way, into an infrared and ul- 
traviolet part. The ultraviolet part is modeled 
by an effective action made up of small Wilson 
loops, while the infrared part is determined pre- 
cisely using a Lanczos procedure. To date, they 
have generated a sequence of gauge fields using 
just the infrared part of the determinant and have 
found that the exact contribution of the ultravi- 
olet part is well fit by an expansion in a small 
number of Wilson loops. A somewhat surprising 
result also comes from the generation of the gauge 
fields: they first update links and then put in a 
global accept /reject step based on the infrared 
part of the determinant. They find reasonable 
acceptance for this, even for large volumes. 

2.3. Multiboson Algorithm 

de Forcrand and collaborators |^ have been 
testing an evolved version of Luscher's multibo- 
son algorithm |^ in a real QCD context and com- 
paring its performance to a standard, but highly 
tuned, hybrid Monte-Carlo code of the SESAM 
collaboration. This new work uses only 24 bo- 
son fields for a simulation with m-n /vip = 0.833, 
a dramatic improvement. This is achieved by an 
"ultraviolet filtering" of the determinant, where 
the determinant is filtered by an exponential term 
made up of weighted small Wilson loops, and a 
cancelling term enters the gauge action. (This 
filtering is similar to that used by 0|, but here 
is part of an exact algorithm.) An optimized 
quasi-heatbath is then applied to the combined 
boson-gauge system. They conclude that for a 
\& X 24 lattice, with (3 = 5.6 and the equivalent 
to K = 0.156 for Wilson fermions, the multiboson 
algorithm decorrelates the plaquette better than 
HMC, measuring both algorithms in units of ap- 
plications of the Dirac operator. 



3. SIMULATION KEY 

To distinguish the different actions, the abbre- 
viations below will be used (following the CP- 
PACS collaboration). 



gauge action 
P 
R 

fermion action 
S 

W 

c 

D 



plaquette 
RG improved 



staggered 
Wilson 
Clover 
domain wall 



The RG improved gauge action is the one pro- 
posed by Iwasaki, et. al. 0. 

In the tables of simulations parameters, three 
results are generally listed: mpO, mj^/mp and 
TTipL, where L is the spatial lattice extent and 
a is the lattice spacing. The value of rupa (where 
available) is for the quark mass which gives the 
quoted rriT^/mp. Also note that a 770 MeV p in 
a 3 fm box has nipL = 11.6. Since these tables 
should only serve as a guide to the general fea- 
tures of the data sets, errors are not listed and 
only 2 significant figures are given. Please see the 
original work for more details. 

The hybrid Monte Carlo (HMC) algorithm has 
been used for the 2 flavor Wilson and domain 
wall and 4 flavor staggered simulations. For 2 
flavor staggered simulations the 'R' algorithm of 
Gottlieb, et. al. ^ has been employed. The tra- 
jectory counts in the tables are from the groups 
themselves and the definition of a trajectory dif- 
fers between the groups. 

4. STAGGERED FERMIONS 

Recent staggered simulations have been done 
by the MILC and Columbia groups and are listed 
in Table |l|. MILC reported hadron spectrum re- 
sults last year p| a nd new results for decay con- 
stants this year |lCf| . They found a continuum ex- 
trapolation of mN/iTip for staggered fermions in 
full QCD gives a larger value than in the quenched 
continuum. The Columbia group has new data 
for 2 and 4 flavor QCD with staggered fermions, 
with simulations still underway pi|. 
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Table 1 

Summary of staggered simulations. 

MILC - PS action - 24^ x 64 - 2 flavors 



0.8 



(3 m 


nipa rUpL 


m^/nip 


traj. 


5.6 0.08 


0.98 23 


0.76 


2000 


5.6 0.06 


0.87 21 


0.74 


2000 


5.6 0.04 


0.75 18 


0.71 


2000 


5.6 0.02 


0.59 14 


0.63 


2000 


5.6 0.01 


0.50 12 


0.53 


2000 


CU QCDSP 


- PS action - 


16^ X 32 - 


2 flavors 


f3 m 


TUpa nipL 


m-^/mp 


traj. 


5.7 0.015 


0.48 7.7 


0.63 


5000 


CU QCDSP 


- PS action - 


16^ X 32 - 


4 flavors 


f3 m 


mpfl rupL 


niTr/mp 


traj. 


5.4 0.02 


0.50 8.0 


0.71 


5180 


5.4 0.015 


0.48 7.7 


0.67 


4750 


CU QCDSP 


- PS action - 


24^ X 32 - 


4 flavors 


[3 m 


mpfl rUpL 


m^/nip 


traj. 


5.4 0.02 


0.49 12 


0.72 


5000 


5.4 0.01 


0.37 8.9 


0.66 


5000 



There are existing results from 2 flavor stag- 
gered simulations at /3 = 5.7 and m = 0.01 which 
give some information about run lengths. Figure 
shows, starting at the top, niN^ nip and from 
the Columbia group ^ and Fukugita, et. al. [|l3|. 
The masses labeled "Ik" are from 1,000 trajecto- 
ries on a 20'* lattice, "3k" from 3,000 trajectories 
on 16^ X 32 and "10k" from 10,000 trajectories 
on 16^ X 40. The points to the left of "10k" are 
the 10,000 trajectory run broken up in to nine 
1,000 trajectory runs (plus thermalization) [Q. 
The solid horizontal lines are the "10k" masses 
and the dashed horizontal lines give ±5% of the 
central value. 

The similar size for the errors on the 3k and 
10k runs (the same analysis was used for both) 
indicates that long correlation times were likely 
missed in the 3,000 trajectory run. The larger 
error bars for the 1,000 trajectory runs may be 
due to short term noise effects. For these masses 
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Figure 1. A comparison of hadron masses for 
a 10,000 trajectory run (10k), a 3,000 trajectory 
run (3k) and a 1,000 trajectory run (Ik). 



and couplings, achieving reliable errors at the few 
percent level likely requires more than 10,000 tra- 
jectory run. 

Another important question for full QCD sim- 
ulations is the volume required. We now investi- 
gate this question through the splitting between 
parity partners in the hadron spectrum. We will 
see that, at least for staggered fermions, this is a 
sensitive indicator of finite volume effects. 

MILC studies of finite volume effects in 2 flavor 
dynamical staggered simulations for couplings up 
to (3 — 5.6 show lattice volumes of 2.5-3 fermi 
remove finite volume effects for m,r/wp ^ 0.5. 
Their most recent (3 — b.& results for a 24* x 
48 lattice |l^] are shown in Figure There is 
no sign of parity doubling, consistent with their 
statement about finite volume effects. 

Fukugita, et. al. also studied finite volume ef- 
fects but at weaker coupling {(3 = 5.7). Their 
data is plotted in Figure ^, revealing parity dou- 
bling in the m ^ limit for the N and N' . 
(The N' is the staggered fermion parity partner 
of N .) Finite volume effects are likely distort- 
ing the baryons, but not the mesons. Which way 
m^v and mjv move for larger volume is an open 
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Figure 2. Some of the MILC 24^ x 48, 2 flavor 
staggered spectrum plotted versus quark mass. 
The lines merely connect the points. 



Figure 3. Staggered 2 flavor simulations on 20'^ 
lattices showing parity doubling in ttiat and iriN' 
as m 0. 



question, currently being addressed by 24"^ simu- 
lations underway at Columbia. 

The Columbia group has previously reported 
parity doubling for 4 flavor QCD on a 16"^ x 32 
volume with /3 = 5.4. (These parameters give a 
rho mass within a few percent of the rho mass for 
2 flavors on a 16'^ x 32 lattice.) Figure^ shows our 
results, where the m = 0.01 point is from the 256- 
node Columbia machine, the m ~ 0.015 point is 
from QCDSP and the m = 0.02 result was calcu- 
lated on both the 256-node machine and QCDSP, 
which agree within errors. Parity doubling is clear 
for both mesons and baryons as m — > 0. 

The Columbia group has now completed a 
5,000 trajectory simulation on a 24"^ x 32 volume 
usmg QCDSP. Figure | shows that the degener- 
acy between the hadrons in the m —^ has gone 
away. There is very little change in the m = 0.02 
masses, but for m = 0.01 nip and mjv have 
dropped by ~ 20% (nip: 0.438(8) ^ 0.373(6), 
tun: 0.690(21) ^ 0.574(9)). This clearly shows 
finite volume effects distorting the parity split- 
tings and here it primarily effects the nucleon and 
rho masses. Also, given the large decrease in nip, 
the larger 24^ lattice has nipL = 8.9 very close to 



7.0 for 16^. 

Four flavors likely make this finite volume effect 
more pronounced, but it is a warning for 3 flavor 
simulations. We are currently simulating with 2 
flavors on a 24'^ x 32 volume to sec how large 
the effects are there. Unfortunately, these large 
finite volume effects are masking any information 
about the role of the determinant in the parity 
splittings. 

A final message about run lengths from the 4 
flavor staggered simulations is shown in Figure 
1^. The upper line is the pion propagator at a 
distance 10 lattice sites from the source for the 
m — 0.01 simulation, plotted against trajectory 
number. The lower line is the same propagator 
for the m = 0.02 simulation. The m — 0.01 simu- 
lation shows fluctuations on a few thousand tra- 
jectory time scale. This is clear evidence that very 
long runs are needed as the quark mass is made 
smaller. (Long autocorrelation times for topolog- 
ical charge for 4 flavor staggered simulations on 
a 16^ X 32 volume with f3 = 5.35 and m = 0.01 
have also been seen 0.) 
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Figure 4. Staggered 4 flavor simulations on 16^ x 
32 lattices showing parity doubling in iriN, uin' 
and TUp, TTia^ as rn ^ 0. 
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Figure 5. Staggered 4 flavor simulations on 2A^ x 
32 lattices showing parity doubling eliminated for 
larger volumes. 



5. WILSON FERMIONS 



The SESAM |TD, UKQCD § and CP-PACS 
|2| collaborations reported on full QCD simula- 
tions last year and UKQCD ^ Q and CP- 
PACS (21) Q have new results this year. The 
run parameters are given in Table |^ and ^ Both 
UKQCD and CP-PACS are using clover improved 
Wilson fermions; UKQCD uses Csw determined 
by the Alpha collaboration and CP-PACS uses a 
tadpole improved value. 

UKQCD chooses (3 and k to keep the lattice 
spacing, as determined using the Sommer scale 
ro and the static quark-anti-quark potential, con- 
stant. The physical volume then corresponds to 
1.7 fm for all lattice spacings they consider. As 
discussed above for staggered fermions, this vol- 
ume can be expected to be rather small. They 
do report evidence that the potential at small r 
for the dynamical simulations lies below the value 
for quenched simulations. In addition, plotting 
vector meson masses against pseudoscalar me- 
son masses, they see a trend toward the physical 
{K,K*) value. 
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Figure 6. The pion propagator at separation 
10 for 4 flavor simulations on 24'^ x 32 lattices 
plotted versus trajectory number. Much longer 
autocorrelation times seem visible for the upper 
plot (rn = 0.01) when compared to the lower 
(m = 0.02). 



6 



Table 2 

Summary of UKQCD and SESAM 2 flavor Wil- 
son fermion simulations. 

SESAM - PW action - /? = 5.6 - 16^ x 32 



K nipa TUpL rriT^/nip traj. 

0l56 (153 8^5 (183 5000 

0.1565 0.50 8.0 0.81 5000 

0.157 0.46 7.4 0.76 5000 

0.1575 0.41 6.6 0.68 5000 



Table 3 

Summary of CP-PACS 2 flavor Wilson fermion 
simulations. 

CP-PACS - RC action - /3 = 1.80 - 12^ x 24 









traj. 


0.1409 


1.60 


0.81 


6250 


0.1430 


1.60 


0.75 


5000 


0.1445 


1.60 


0.70 


7000 


0.1464 


1.60 


0.55 


5250 



UKQCD - PC action - 16^ x 32 
13 = 5.29, 5.26, 5.2, 5.2, respectively 



CP-PACS - RC action - /3 = 1.95 - 16^ x 32 



K 




rUpQ 


TTlpL 


m^/nip 


conf. 


0.1340 


1.92 


0.70 


11 


0.83 


101 


0.1345 


1.95 


0.65 


10 


0.78 


101 


0.1350 


2.02 


0.59 


9.4 


0.69 


150 


0.1355 


2.02 






0.58 


102 



The CP-PACS collaboration has been involved 
in extensive simulations of full QCD with a va- 
riety of lattice spacings, quark masses and vol- 
umes. Their parameter choices keep the spatial 
size fixed at ~ 2.4 fm, using the p mass at the 
physical rriT^/mp value to set the scale. With their 
full QCD data set, they can extrapolate to the 
continuum, with fixed finite volume effects. To 
date, they only have 2,000 trajectories for their 
/3 = 2.10 point, which experience with staggered 
fermions suggests may not be long enough. They 
are addressing this issue. 

CP-PACS found that a major difficulty with 
the quenched hadron spectrum is the failure of 
a unique strange quark mass to give the physi- 
cal values for the K and masses. They have 
addressed this question with their new data and 
one of their results is shown in Figure 0. They 
find evidence that the a — » extrapolation of 
the lattice result is much closer to the physical K 
mass (o) than the quenched a — > value. Extrap- 
olations of the octet and decuplet baryon masses 
[|2l[| have larger errors and definitive conclusions 
cannot be drawn. 

Given their evidence that a single value of the 
strange quark mass determines both the K and 
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7000 
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7000 
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0.69 




7000 
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7000 


CP-PACS - RC action - 


/3 


= 2.10 - 243 X 48 
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niTr/mp 




traj. 


0.1357 
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2000 


0.1367 
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2000 
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2000 
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(f) masses, they have also determined the a ^ 
strange quark mass at ^ — 2GeV. One of their re- 
sults is shown in Figure ^. They find various ways 
of determining the strange quark mass all agree 
in the a — > limit and there is a large difference 
between the quenched and unquenched values for 
this mass. Their final result is rus = 87(11) MeV 
(0 input) and nis = 84(7) MeV {K input). These 
values are considerably lower than those from 
phenomenology. They also find mu,d = 3.3(4) 
MeV. 

The CP-PACS collaboration also reported a 
value for the flavor singlet pseudo-scalar meson 
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Figure 7. CP-PACS meson masses using the 
mass to set the strange quark mass. 
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for 2 flavor QCD, referred to as the rj [g2|. Us- 
ing a volume source without gauge fixing (the 
Kuramashi method) to measure the disconnected 
quark diagrams, the mass difference between the 
TT and rj can be calculated, at least when the 
disconnected diagrams are not far separated in 
time. They find m,, non-zero in the chiral limit 
and a subsequent continuum extrapolation gives 
TO,, = 863(86) MeV. 

UKQCD has also reported a preliminary value 
for the 77 mass for their set of dynamical lat- 
tices. They employ a variance reduction tech- 
nique in finding the propagator for the discon- 
nected diagrams and get an r) mass around 800 
MeV in the chiral limit, with an uncontrolled sys- 
tematic error. 

Both groups see little problem with autocor- 
relation times for these topologically sensitive 
measurements. UKQCD finds an autocorrelation 
time smaller than that found by the SESAM col- 
laboration. 

6. DOMAIN WALL FERMIONS 

The development of domain wall and overlap 
|23| fE4l [G5| formulations of lattice fermions has 



produced a great deal of excitement for both an- 
alytical and numerical studies. These formula- 
tions offer a way of separating the chiral limit 
from the continuum limit and may lead to a lat- 
tice formulation of non-abelian chiral gauge theo- 
ries For QCD-like theories, quenched simu- 
lations with domain wall fermions were first done 
by and were reviewed last year j28j . Concur- 
rently, simulations of the Schwinger model [ p9[ , 
including dynamical fermions, showed that do- 
main wall fermions produced the expected physics 
and were compatible with standard HMC algo- 
rithms. 

Two areas where using fermions with better 
chiral properties are of particular interest are in 
simulations studying QCD thermodynamics |30) 
and matrix elements 1 31 1 |^ |Q . The character 
of the QCD phase transition is governed by the 
chiral symmetries of the theory. For matrix ele- 
ment calculations, chiral symmetry can be vital 
for controlling operator mixing and using chiral 
perturbation theory as a guide. Both of these 
areas are under active study with the QCDSP 
computers, using domain wall fermions. 

The Columbia group has been studying 
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full QCD thermodynamics with domain wall 
fermions. The studies to locate the transition re- 
gion and set the parameters to use are detailed in 
. As part of this work, zero temperature scale 
setting calculations are also required, which we 
will focus on here |^^. The full QCD, zero tem- 
perature domain wall simulations done to date are 
detailed in Table ||. For these dynamical simula- 
tions, a heavy bosonic field, frequently called the 
Pauli-Villars field, is needed to remove the bulk 
infinity that occurs when the extent of the fifth 
dimension is sent to infinity. 

Table 4 

A summary of 2 fiavor domain wall fermion simu- 
lations. All simulations use a domain wall height 
of 1.9. 

CU QCDSP - PD action - (3 = 5.325 - 8^ x 32 
m Ls mpQ rUpL mj^/mp traj. 
~0M 24 TS TOA 064 1170 
0.02 24 1.2 9.6 0.55 955 

CU QCDSP - PD action - f3 ^ 5.325 - 16^ x 16 
m Ls mripa rUpL mT^/mp traj. 
~0m 24 El 06 057 560 

CU QCDSP - RD action - /3 = 1.9 - 8^ x 32 

m Ls rripa rUpL rriT^/mp traj. 
~O02 24 L2 06 052 875 

CU QCDSP - RD action - = 2.0 - 8^ x 32 
m Ls nipa rUpL mT^/mp traj. 
~O06 24 O 8^8 066 960 
0.02 24 0.95 7.6 0.50 1010 
0.02 48 1.0 8.0 0.43 760 



Scale setting calculations to support thermo- 
dynamics studies are necessarily at fairly strong 
coupling. A first positive result is that even for 
these course lattices, the HMC algorithm exhibits 
no problems thermalizing and evolving lattices. 
A clear region where the Wilson line and the chi- 



ral condensate undergo a rapid crossover Q is 
seen and scale setting calculations have been done 
there. An important question is how the resid- 
ual quark mass, mics (due to mixing of the light 
modes between the two walls), depends on the 
extent of the fifth dimension. 

For the PD action, the critical coupling on 
an TVt 4 lattice is ^ = 5.325, for Ls = 24, 
ruf = 0.02 and M = 1.9. Figure || shows hadron 
masses for this /? determined on 8^^ x 32 lat- 
tices for dynamical fermion masses of 0.02 and 
0.06. By extrapolating to zero, one finds 
rrires = 0.059(2). (A similar value, within errors, 
can also be extracted from the Ward-Takahashi 
identity |^.) This is clearly a large mass, com- 
pared to the input mass of 0.02. 

Extrapolating the rho and nucleon mass to the 
point where vanishes gives — 1.02(7), 
niN = 1.14(9) and rUp/mN = 1.10(12). (The 
errors are calculated by simple propagation of 
errors.) While the length of the simulations 
(~ 1,000 trajectories) may be suflacient at these 
strong couplings, only a single volume and lattice 
spacing has been studied. However, it is encour- 
aging that this ratio is much closer to its physical 
value than for other fermion methods at this lat- 
tice spacing. 

It is important to study the effects of increasing 
the fifth dimension on the residual mass. Figure 
[To| shows the residual mass versus Ls (the length 
of the fifth dimension) for simulations on 8"^ x 4 
lattices with the PD action at /? = 5.2 with a 
quark mass of 0.02 and a domain wall height of 
1.9. These lattices are in the confined phase. The 
residual mass is determined from the axial ward 
identity, as discussed in . The function shown 

ismj.™°^^ = 0.17(2) xexp(-0.026(6)xL3)- The 
data shows a vanishing residual mass in the large 
Ls limit, but the falloff is very slow. Increasing 
the fifth dimension by 24 drops the residual mass 
by about a factor of 2. 

In an effort to decrease mres for a fixed Ls, 
the Columbia group has studied the renormaliza- 
tion group improved action proposed by Iwasaki. 
For quenched simulations, Wrcs is made smaller 
for a given Ls, but the behavior with Ls may 
not be improved iQ. Results for hadron masses 
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Figure 9. Columbia QCDSP hadron masses for 
simulations with dynamical domain wall fermions 
with the PD action. The scale of the phase tran- 
sition on an A^t = 4 lattice determines the param- 
eters. 
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Figure 10. Wres from the Ward-Takahashi iden- 
tity as measured on 8^ x 4 lattices in the confined 
phase with the PD action and /3 — 5.2. 



for dynamical simulations with the RD action are 
shown in Figure ^ These were carried out on 
8^ X 32 lattices for /? = 2.0 with = 24. From 
the behavior of ra^ one finds TOrcs = 0.013(2). 
At the point where vanishes rup = 0.855(49), 
ruN = 1.07(14) and mp/mN = 1.25(14). (Once 
again, these are naive errors.) 

While mrcs is smaller for the RD simulations at 
/3 = 2.0 than the PD simulations aX (i ~ 5.325, 
the physical lattice scales are different. For the 
RD action, the Nt = A phase transition \s aX (3 — 
1.9. We do not have two dynamical masses for the 
RD action at /3 = 1.9, but we can compare the 
pion masses at m = 0.02. For the PD action at 
(3 = 5.325, = 0.654(3) and for the RD action 
at /3 = 1.9, TOjr = 0.604(3), a slight improvement. 
(However, the uncertainty in the determination 
of the critical temperature could easily account 
for this difference.) 

A direct comparison of the pion mass for the 
RD action (/3 = 2.0, 8^ x 32) for Ls = 24 
gives TOtt = 0.475(7), while for Ls = 48, m^r — 
0.420(10). mros is decreasing as Ls is increased, 
although the rate is slow. 



7. CONCLUSIONS 

With the current- Tcraflops scale computers, 
full QCD simulations have reached the point 
where the finite volume, long simulation time, 
small quark mass and continuum limits can be 
probed, but not concurrently. The HMC and 
HMD algorithms continue to be the techniques 
of choice, but the multiboson algorithm has been 
evolved to a competitive level. Staggered 4 flavor 
QCD shows large finite volume effects, even for 
m^/mp > 0.5 and the parity splittings for light 
hadrons are quite volume sensitive. For currently 
accessible weaker coupling simulations, runs of 
length greater than 10,000 trajectories are prob- 
ably needed. 

CP-PACS has improved their results on the 
continuum limit of 2 flavor QCD and find contin- 
uum meson masses closer to the physical values 
than for the quenched case. This leads them to 
a determination of the strange quark mass which 
is lower than that expected phenomenologically. 
For staggered fermions, the new 2 flavor points 
being produced by Columbia will augment the 
existing continuum extrapolation of the MILC 
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Figure 11. Columbia QCDSP hadron masses for 
dynamical domain wall simulations with the RD 
action for /? = 2.0. 



group. 

Full QCD domain wall fcrmion simulations are 
straightforward, but require a large fifth dimen- 
sion at strong coupling. It is encouraging that on 
very coarse lattices tojv/top is much lower than 
for other fcrmion formulations. Further studies 
will be needed to check the scaling properties for 
domain wall fermions to see if the use of coarser 
lattices can offset the extra calculational cost of 
the fifth dimension. 
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